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Abstract. We calculate the wake induced in a hot QCD plasma by a fast parton in 
the framework of linear response theory. We discuss two scenarios: (;) a weakly 
coupled quark-gluon plasma described by hard-thermal loop perturbation theory and 
(//) a strongly coupled quark-gluon plasma which resembles a quantum liquid. We 
show that a Mach cone can appear in the second scenario, but not in the first one. 



1. Introduction 

The Oxford English Dictionary records three different meanings of the noun wake: 

• a trail of disturbed water or air left by the passage of a ship or aircraft; 

• a watch or vigil held beside the body of someone who has died; 

• (especially in Ireland) a party held after a funeral. 

Since the Greek word aX)jU7tocnov describes a gathering of friends to eat, drink, and con- 
verse, the Symposium in memory of our deceased friend and colleague Gerhard Soff cer- 
tainly fits the third definition. My lecture will be mainly concerned with the first definition, 
in a slightly generalized sense. To define the appropriate context, we first need to review 
some of the salient results from the experiments conducted at the Relativistic Heavy Ion 
Collider (RHIC). 

The quenching of QCD jets in relativistic heavy ion collisions due to the energy loss 
suffered by hard partons as they traverse dense matter was proposed more than two decades 
ago as an important indicator for the creation of a quark-gluon plasma [[D El |3. Over the 
past five years, this phenomenenon has been extensively studied experimentally at RHIC. 
One looks for hadrons with a high transverse momentum pT, which are produced when 
an energetic, hard scattered parton fragments into hadrons far outside the remnants of the 
nuclear collision. Since both, the scattering probability and the fragmentation probability, 
are known or calculable in QCD, the sole unknown is the amount of energy lost by the 
parton before it fragments into hadrons. This makes it possible to deduce the energy loss 
from the measured hadron yield. 

The RHIC data clearly demonstrate a strong suppression of the emission of high-pj- 
hadrons in Au + Au collisions. The suppression is found to grow in severity with increasing 
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centrality. The suppression factor Raa, defined as the ratio of the hadron yield in Au + 
Au collisions compared to the yield in p + p collisions scaled by the appropriate number 
of binary nucleon-nucleon interactions in the nuclear collision, reaches about 1/5 in the 
most central events [|4l. Most of the observed hadrons originate near the surface of the 
interaction region oriented toward the detector. Hadrons from the companion jet emitted in 
the opposite direction, which has a much longer path through the medium, are even more 
strongly suppressed [jsj. 

The main emphasis in theoretical studies of jet quenching has been on the description 
of the energy loss which the leading parton suffers due to the emission of a secondary 
partonic shower when traversing the medium. Reviews of the theory of medium induced 
energy loss and its associated phenomenology are found in [T 'W, "5" '10]. The main result 
is that the energy loss of a hard parton in dense QCD matter is dominated by radiative 
processes involving gluon emission after elastic collisions of the parton with color charges 
(mostly gluons) contained in the medium. Coherence effects lead to a quadratic dependence 
of the total energy loss on the traversed path length L. The stopping power of the medium 
is encoded in the quantity 



where p denotes the gluon density of the medium and da/dq^ is the differential cross 
section for elastic scattering. 

An interesting question is: What happens to the radiated energy? Or more generally: 
How does the medium respond to the penetration by the hard parton? The RHIC experi- 
ments are just beginning to address this intriguing question. First indications of a medium 
response have been seen in three sets of data: 

1. Although energetic hadrons in the direction opposite to a high-y^j- trigger hadron are 
almost completely suppressed in central colhsions [|5j|, one finds an increased yield 
of soft, low-pT hadrons [[H] 1121 . 

2. The angular distribution of the soft hadrons emitted in the direction opposite to an 
energetic hadron does not seem to peak at 180° but at a smaller angle [[21 1141 . 

3. Correlated two-hadron emission in the same direction in the momentum range 0.15 
GeV/c < pt <2 GeV/c is enhanced [ 15l. 

Since the quark-gluon plasma is a plasma, after all, it is natural to begin a theoretical 
investigation of the medium response by applying the tools which have been successfully 
used to describe the response of a metallic electron plasma to the penetration by a fast 
ion. Such a charged projectile induces a wake of charge and current density in the target, 
accompanied by induced electric and magnetic fields. The wake, which has the shape of a 
Mach cone, reflects significant aspects of the response of the medium. 

After a brief review of the evidence for this phenomenon in condensed matter physics, 
we apply the same methods of linear response theory to the system of a relativistic color 
charge traveling through a QCD plasma. We calculate the plasma reponse to an external 
point source traveling at a velocity close to the speed of light. In this framework, quantum 
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effects are included implicitly via the dielectric functions, and £7-. For simplicity, we 
shall assume that the medium is homogeneous and isotropic, and we disregard finite size 
effects. We then consider two models for the color response of the QCD plasma: (/) the 
response predicted by QCD perturbation theory, and (//) a response of the kind expected for 
a strongly coupled, "liquid" plasma. We present the results of our calculations of the wake 
structure behind a fast color charge for both scenarios and discuss the general conditions 
which must be met for the wake to have a Mach cone-like structure [[l6| 1171 . 



2. Plasma linear response theory 



The linear response of the plasma to an external electromagnetic field has been extensively 
studied in plasma physics (see e. g. ['181). In this formalism, a dielectric medium is charac- 
terized by the components of the dielectric tensor e,y(a),A:). For an isotropic, homogenous 
medium the dielectric tensor can be decomposed into its longitudinal and transverse com- 
ponents characterized by the dielectric functions ei{(x),k) and e7-(co,fe); 



(2) 



Here V^jj = kikj/k^ and fFj- = 1 — fP/, are the longitudinal and transverse orthonormal pro- 
jectors with respect to the momentum vector k. One can relate the dielectric functions Ez, 
and £7 to the self-energies and of the in-medium photon via [ .19J : 



nL{&,k) 



£T{(0,k) = 1 - 



UT{(a,k) 



(3) 



Using Maxwell's equations and the continuity equation in momentum space, the total elec- 
tric field £tot in the plasma is related to the external current y'ext via: 



k^ 



4%^ 

Etot{&,k) ^ — ;ext(«,^) 

103 



(4) 



Equation @ has propagating solutions when the determinant constructed from the elements 
of the tensor vanishes: 

k^ 



det 



0. 



(5) 



This equation governs the dispersion relation for the waves in the medium. It can be diag- 
onalized into purely longitudinal and transverse parts, yielding dispersion relations for the 
longitudinal and transverse dielectric functions [ 1 8 1 : 



et(co,^) =0, 



ET{03,k) - {k/03f 



(6) 



These equations determine the longitudinal and transverse plasma modes. The longitudinal 
equation is the dispersion relation for density fluctuations in the plasma, namely space- 
charge fields which can propagate through the plasma without far away from an external 
perturbation. 
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The charge density induced in the wake by the external charge distribution is: 

Pind = (^-^ - 1^ Pext- (7) 

In the transverse gauge, the induced charge density is related to the induced Coulomb 
potential via the Poisson equation: k^(^ind = 47tpind- Since one can relate the total electric 
field to the induced charge in linear response theory by 

Jind = J^(l-e)£tot, (8) 

471 



a direct relation between the external and the induced current can be derived using 0: 

/ext. (9) 



7ind = 



1 \ 0)2(1 -Er) 



The induced charge and the induced current obey the continuity equation: 

ik ■ jmd - /COPind = 0. (10) 

Finally, we need to specify the form of the external current. For a fully stripped ion 
(or a single energetic parton) it is appropriate to assume the current and charge density of 
a point-like charge moving along a straight line trajectory with constant velocity v, whose 
Fourier transform is given by 1 1201 : 

/ext 2%qv?>{(a-v-k), 

Pext = 27t^5(CO-V-^). (11) 

All equations given above, from Q to M \l . immediately generalize to QCD by the 
simple addition of a color index a = 1 , . . . , 8 to the charge, the current, and the field strength, 
which are all in the adjoint representation of color SU(3) : 

q^q", (PJ)-(P",/'), E^E". (12) 

Since we are limiting the treatment of the medium response to effects that are linear in the 
perturbation, all nonlinear terms arising from the nonabelian nature of the color field are 
discarded. The strength of the color charge of the projectile is defined by q^'q" — C2OC4 with 
the strong coupling constant a, ~ g'^/4-Ti and the quadratic Casimir invariant C2 (Cp =4/3 
for a quark or antiquark and d = 3 for a gluon). In this linearized treatment one disregards 
changes of the color charge while the particle is propagating through the medium. The non- 
abelian character of the QCD plasma only enters indirectly via the chromo-dielectric func- 
tions ei/j{(x),k), which are scalars in color space [|^ |22l |23 1- Due to the self-interacting 
nature of the gluon in SU(3), gluons contribute to the polarization of the medium; in fact, 
they make the largest contribution. 
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The non-radiative part of the energy loss of the incident (color) charge is given by the 
back reaction of the induced (chromo-)electric field onto the incident particle. The energy 
loss per unit length is given by [ 1 8 1 : 



— = q-KaEini{x = vt,t), 
dx V 



(13) 



where the induced electric field Ei^d is the total electric field minus the vacuum contribution. 
Using the inverse of @, the induced field is given by: 



^ind ■ 



1 

El 



1 hPL 



CO"- 



4% 
m 



Je\t- 



From jl3> and (I14> the non-radiative energy loss per unit length is given by [ 1231 : 



dE 

dx 



a\-l^ 



(14) 



(15) 



where (sik = v -k. 



3. Quick flashback: Electron wakes in metal foils 

The formation of electromagnetic wakes induced by a fast ion in the electron plasma of a 
thin metal foil was investigated starting around 1980 by Groeneveld and collaborators at 
Frankfurt [ 24 25 1. Their study was motivated, in part, by calculations [26 27) done by a 
graduate student, Wolfgang Schafer, ' who had solved the equations of the previous section 
for a dielectric function of the Bloch type [ 1281 

s. = l + ^^^ {k<k.). (16) 

where cOp is the plasma frequency and u denotes the sound velocity in the plasma. The 
calculations showed that the wake in the electron plasma has the form of a Mach cone, 
similar to the Mach shock phenomena predicted for relativistic nucleus-nucleus collisions 
r i29l 1301 . In our publications, we not only showed the detailed shape of the spatial distri- 
bution of induced plasma charge and current (see Fig.[Q, but we also predicted the angular 
distribution of the electrons which would be emitted when the current wave hits the surface 
of the metal foil (see Fig[3. The peak of this distribution roughly coincides with the "Mach 
angle", given by 

cpM = arccos(M/v). (17) 

The first published data confirmed the presence of directed electron emission, and the 
dependence of the peak angle on the beam energy qualitatively agreed with our predictions 
r i24l . However, it turned out that the results were not well reproducible unless the surface of 

'Before doing this work for liis doctoral tliesis, Wolfgang had worked under Gerhai'd Soff's supervision cal- 
culating the influence of the vacuum polarization potential on the motion of two colliding heavy nuclei. 
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Fig. 1. Left: Spatial distribution of the induced charge density behind an ion traveling 
through a metal. The plot represents a radial cut through the Mach cone. Right; Predicted 
angular distribution of the electrons contained in the wake current. 



the foil was extremely clean. Groeneveld's team therefore later repeated the measurement 
with foils whose surface had been sputter-cleaned. The new data (see Fig.|2ll showed much 
more pronounced and reproducible peaks [ '25V Their positions as a function of beam 
energy and Fermi energy of the target nicely followed the Mach relation (I17> in agreement 
with the predictions \ \27l . 

The experimentalists continued to study many aspects of this phenomenon in great 
detail, including the energy spectra and angular distributions of the ejected electrons, the 
refraction of the Mach wave at the planar surface of the foil, and the response of high Tc 
superconductors ['311. In recent years, the collective plasma waves excited by fast ions in 
soft biological tissue have been considered as a mechanism that contributes to the damage 
to living cells by fast ions [ 32|, which is an important factor in cancer therapy with 
heavy ion beams [j33J. 



4. Color wake in the high temperature approximation 

We now return to the problem of interest to us: the medium response to an energetic parton. 
We shall discuss two qualitatively different scenarios. We first assume that the plasma 
is in the weakly coupled high temperature regime, where the gluon self-energy can be 
described by the leading order of the high temperature expansion, T ^ (0,k, commonly 
called the hard-thermal loop (HTL) approximation [|34] 1351 . This regime can be expected 
to be realized far above the deconfinement temperature Tc. The dielectric functions read [ 
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Fig. 2. Left; Angular distribution of electrons emitted after bombardment of thin metal foils 
by fully stripped ions. Right: Dependence of the peak position in the angular distribution 
for various target materials (C, Al, Cu) and different beam energies, in comparison with the 
predicted Mach angle for a carbon target (from r.25J). 



where x — &/k. We first discuss the induced charge and current densities and yind- The 
Fourier transform of Q in cylindrical coordinates is given by: 



mlq" 



{2%fv Jo 
d(i) exp 



t/K'K'7o(Kp) X 



(20) 



m 



1 

El 



where k — -^/k^ + a)^/v'2 and CO = m^co', k = m^K'. This shows that the induced charge 
density p"^^ is proportional to m^. The cylindrical symmetry around the jet axis restricts 
the form of the current density vector j"^^. It only has non-vanishing components parallel 
to the beam axis: 



23 



(27r) 



(iK'K'7o(Kp) 



lico'exp 
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1 -Er 
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and radially perpendicular to the beam axis: 



dydydJx (Kp) 



c/co'exp m{^~t^ 
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Fig. 3. Spatial distribution of the induced charge density for a color charge traveling with 
velocity v/c = 0.99 in a high temperature QCD plasma where the HTL approximation 
applies. The left plot shows equi-charge density lines in the rest frame of the charge. 

Again, the components of the current density are proportional to m^,. 

For the dispersion relations ( I18> . longitudinal and transverse plasma modes can only 
appear in the time-like sector of the co,^ plane [[19] 1211 . Therefore, collective excitations 
do not contribute to the charge and current density profile of the wake. Emission analogous 
to Cherenkov radiation and Mach cones are absent, but the charge carries a screening color 
cloud along with it. Fig.|3] which shows the charge density of a colored parton traveling 
with V = 0.99c, illustrates this physically intuitive result. 

In spite of the absence of a Mach cone, the particle loses energy due to elastic colli- 
sions in the medium, which can be described by formula (I15> . This mechanism of energy 
dissipation has been studied in [ i23J . The integrand in M5\ contributes to the integral in the 
space-like region only, where < 1, and hence does not get contributions from frequen- 
cies where collective plasma modes exist. This is consistent with the fact that such modes 
are not excited by the moving color charge. 

5. Charge wake induced in a strongly coupled QGP 

In the second scenario we investigate what happens if the QCD plasma is a strongly coupled 
quark-gluon plasma (sQGP) that can be described as a quantum liquid. Our consideration 
of the second scenario is motivated by the RHIC experimental results on collective flow, 
which have led to the conclusion that the QCD plasma behaves like a nearly ideal fluid with 
very low viscosity (see e. g. r i36I V This implies that long wavelength collective modes 
are almost undamped, while the short-distance dynamics are strongly dissipative due to 
large transport cross sections. In this respect, the sQGP resembles the electron plasma in 
metallic targets described in the previous section. The idea that a fast parton could excite 
Mach waves in such a QCD plasma and that the sound velocity of the expanding plasma 
could be determined from the emission pattern of the secondary particles traveling at an 
angle with respect to the jet axis was first suggested by Stocker r i37l . The recent work by 
Casalderrey-Solana and Shuryak [ 38 1 explores the formation of a conical flow by a hard 
parton in a quark-gluon plasma within a hydrodynamical framework, but does not specify 
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Fig. 4. Dispersion relation for the plasmon mode for a weakly coupled plasma described 
by the HTL formalism (solid red line) and for a strongly coupled plasma described by the 
modified Bloch formula i24\ . The latter extends into the space-like region (k > co) of the 
(0 — k plane. The dispersion curve of transverse HTL-plasma modes is shown in blue and 
the light cone (co — k) is represented in black. The black dot indicates the point ks where 
the dispersion curve intersects the light cone. 

how the energy of the quenched jet is deposited into the medium. 

There are presently no theoretical methods available for first principle calculations 
of the color response functions in a strong coupled QGP plasma. We therefore confine our 
investigation to a simple model, which encodes the essential differences between a quantum 
liquid scenario and the weak coupling scenario of the last section. The most prominent 
difference is the possibility that a plasmon mode may extend into the space-like region 
of the CO — plane above some threshold value ks- As we already emphasized, the sQGP 
paradigm suggests very low dissipation at small k, but large dissipation at high k. Our 
assumption is that a critical momentum k^ separates the regimes of collective and single 
particle excitation modes in the quantum liquid where the dominant colored modes below 
kc are plasmon excitations with negligible dissipative single-particle coupling. Since we 
are here predominantly interested in collective effects in the plasma, we restrict our study 
to the region k < kc and simply cut off all Fourier integrals at kc- 

To be specific, we assume that the dielectric function of the strongly coupled plasma 
in the k < kc regime leads to a longitudinal dispersion relation of the form: 



where cOp denotes the plasma frequency and m < c is the speed of plasmon propagation, 
here assumed to be constant. In accordance with (123 > we posit the following dielectric 
function 




(23) 
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Fig. 5. Spatial distribution of a induced charge density around a color charge traveling with 
velocity v = 0.55c < u. The left plot shows equi-charge density lines. The density profile 
is similar to that of the cloud surrounding a color charge at rest. 

which differs from the classical, hydrodynamical dielectric function of Bloch [ 28 1 by re- 
maining regular in the limit A:, CO ^ 0. The Bloch function is singular at small k, (O due to 
the mixing of the plasmon mode with the phonon mode. Such a mixing cannot occur in the 
QCD plasma, because the plasmon and phonon belong to different irreducible representa- 
tions (octet versus singlet) of color SU(3) and because the medium is charge symmetric. 

The dielectric function i24i is constructed in such a way that it allows us to study one 
specific aspect of a quantum liquid scenario: that the plasmon mode may extend into the 
space-like region of the 03 — k plane. This behavior is illustrated in Fig.|4]by the dashed 
red line, which cuts through the light cone at kg (at the black dot in Fig.|3 and continues 
into the space-like region. This contrasts with the dispersion relation of the longitudinal 
plasmon mode in the HTL formalism (solid red line), which always remains in the time- 
like region a)(fe) > k. The induced wake structure for a supersonically (in the sense v > u) 
traveling color source in such a quantum liquid scenario is, quite generally, a conical Mach 
wave structure. The principal findings of our study can be expected to hold for any quantum 
liquid with a plasmon branch similiar to i23\ . independent of the exact form of the dielectric 
function. We here assume a speed of plasmon propagation u/c — l/\/3. This differs 
from the speed of plasmon propagation in the small k limit of the HTL approximation, 
u/c ^ \/3/5 \ \2ll . and is also different from the sound velocity in a hadronic resonance 
gasM/c« V0^['39','40l. 

Determining the plasmon mode via Eq. i23\ reveals that the mode is in the space-like 
region of the CO — A: plane for k > k^ = (i>p/ Vc^ — u^ ■ Recall that this is different from the 
high-temperature plasma, where longitudinal and transverse plasma modes only appear in 
the time-like region, |x| — \03/k\ > 1. In the quantum liquid scenario one can expect that 
the modes with low phase velocity \x\ < u/c suffer severe Landau damping because they 
accelerate the slower moving charges and decelerate those moving faster than the wave. A 
charge moving with a velocity lower than the speed of plasmon propagation can only excite 
those modes and not the modes with intermediate phase velocities u/c < \x\ < 1, which are 
undamped [TS* Tl]. In this case, the qualitative properties of the color wake are analogous 
to those of the high temperature plasma: the charge carries a localized screening color cloud 
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Fig. 6. (a) Spatial distribution of the induced charge density from a jet with high momentum 
and fixed color charge q'^ that is traveling with v = 0.99c > u = -y/l/Sc. (b) Plot showing 
equi-charge lines in the density distribution for the situation in (a). 

with it and Cherenkov emission and Mach cones are absent. Using (I20> and restricting the 
integration area to the region k < kc = 2(0p, one can illustrate this for a colored particle 
traveling with v = 0.55c < u. Figure |5] shows the charge density cloud traveling with the 
colored parton. 

If the colored parton travels with a velocity v that is higher than the speed of plasmon 
propagation u, modes with an intermediate phase velocity u/c < \x\ < 1 can be excited. 
The emission of these plasma oscillations induced by supersonically traveling particles is 
analogous to Cherenkov radiation, but different in that the density waves are longitudinal, 
not transverse, excitations of the color field. Figure |6l for a color charge traveling with 
v/c — 0.99 clearly exhibits the emergence of Mach cones in the induced charge density 
with an opening angle given by the Mach relation (I17> . 

We emphasize that the existence of Mach cones is expected in a plasma in general if 
the particle is moving faster than the speed of sound in the plasma and if the dispersion 
relation of the collective mode extends into the space-like region. The wake induced by a 
colored jet in such a setting leads to regions of enhanced and depleted charge density in the 
wake, which have the shape of Mach waves trailing the projectile. 

6. Summary and Outlook 

We have calculated the properties of the color charge density wake of a hard parton travel- 
ing through a quark-gluon plasma in linear response theory for two different scenarios: a 
weakly coupled QGP at T ^ Tc described in the HTL approximation, and a strongly cou- 
pled QGP with the properties of a quantum liquid. We found that the wake in the weakly 
coupled plasma always takes the form of a screening cloud traveling with the particle, while 
the wake in the strongly coupled plasma assumes the form of a Mach cone, if the parton's 
velocity exceeds the speed of plasmon propagation and the collective plasma mode has a 
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dispersion relation extending into the space-like region. 

In general, secondary particle distributions can be used to probe the collective ex- 
citations QCD plasma. If these take the form of Mach cones they should reveal them- 
selves by the directed emission of secondary particles from the plasma, similar to what 
was observed with metallic targets under bombardment by swift ions. If the scenario of a 
strongly coupled QCD plasma is realized in relativistic heavy ion collisions, one can ex- 
pect to observe these cones in the angular distribution of secondary particles associated 
with jets 1381 . As we already mentioned in the introduction, preliminary data from 
the RHIC experiments (see e. g. Fig. 1 in [[2] O) suggest that the angular distribution 
of secondary hadrons in the direction opposite to an energetic hadron is peaked at an az- 
imuthal angle A(|) ^ %. This is different from what is observed in p + p collisions, where 
the maximum associated with the away-side jet is clearly located at A(|) = 7t. In contrast, 
two maxima are located at A(|) w 7i ± 1 . 1 in Au + Au collisions. This phenomenon suggests 
the presence of a Mach shock front traveling with the quenched away-side jet at an angle 
A(|) — 7t ± cpM — 7t ± arccos(M/v). 

If these speculations are confirmed, it will be an interesting theoretical problem to de- 
termine the mechanism that excites the shock front. Possible candidates are the collective 
plasma excitations discussed here [ 16 37 1, the "thunder" effect following localized di- 
rected energy deposition by the quenched jet [ 1381 . or simply the knock-on process due to 
elastic collisions of plasma particles with the away-side hard parton r i41l . 
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